Quantum size effects in thermodynamic superconducting properties of ultra thin films 
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By using multi-bands BCS theory, we have calculated the superconductivity energy gap and the 
critical temperature of a thin-film metallic superconductor. The thermodynamic superconducting 
characteristics such as critical magnetic field, specific heat, as well as the tunneling conductance are 
investigated for varying film thickness and temperature. We find the oscillation of thermodynamic 
superconducting properties with the film thickness, including the thermodynamic critical field H c , 
the specific heat of normal, superconducting state, and position of the differential tunneling conduc- 
tance peak. The two universal constants the nth sub-band energy gap A n at temperature T = OK 
over ksT c , and the specific heat jump at T c over normal state specific heat at T c are independent 
of the film thickness. Their values are the same as in the bulk superconductor. 
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Quantum size effects (QSE) in superconductors have 
been a very attractive topic after some theoretical in- 
vestigations of the superconducting transition tempera- 
ture and energy gap for thin films. [3 Some theoretical 
|H |j] and experimental 0, 0, 0, works have re- 
ported that several quantities were modulated by QSE, 
including electronic structure, electron-phonon interac- 
tion, resistivity, Hall conductivity, work function, surface 
energy and superconductivity critical temperature. In 
a metal film on a semiconductor substrate, the conduc- 
tion electrons are confined by the vacuum on one side 
and the metal-semiconductor interface on the others. 
QSE and the simple stability of metal thin films on 
a supporting; substrate have been discussed in several 
papers. [U El 

The report by Guo et al. is the first definitive and 
quantitative demonstration of T c oscillation with Pb film 
thickness as well as the normal state conductance. E3 Us- 
ing atomically uniform film of lead with exactly known 
numbers of atomic layers deposited on a silicon (111) 
surface, Guo et al. observed oscillations in T c that cor- 
related well with the confined electronic structure and 
recently they demonstrated that film thickness can in- 
deed affect superconductivity behaviors. It has been 
found that when the thickness of a film is reduced to 
the nanometer scale, the film's surface and interface will 
confine the motion of the electrons, leading to the for- 
mation of discrete electronic states known as quantum 
well states (QWS).|lj| QWS have already been observed 
in thin metallic films. [lflj And they change the overall 
electronic structure of the thin film. At a tiny thickness, 
physical properties are thus expected to vary dramati- 
cally with the thickness. Very recently, Daejin et al. us- 
ing a scanning tunneling spectroscopy study to show that 
both energy gap and transition temperature exhibit per- 
sistent oscillation without any suppression at ultra thin 
Pb films (5-18ML).0 



Quantum oscillations can be understood by consider- 
ing QSE in these systems. [l3| The period is fixed for each 
system and equals one half of the bulk Fermi wavelength 
which is related to the average electron density and crys- 
tal structure. E3 Although lots physical properties mod- 
ulated by QSE in thin films have been revealed, the ther- 
modynamic superconducting characteristics such as the 
specific heat, the thermodynamic critical field, tunneling 
conductance as the function of thickness and other pa- 
rameters have not been reported. In this paper, by using 
the multi-band superconductivity theory which improves 
previous theory about QSE in superconductor film, we 
present the results of thermodynamic superconducting 
properties of ultra-thin metallic films. 

From theoretical point of view,[lj the resonance and 
the strong thickness dependence of T c are the character- 
istic features of a thin film superconductivity. The quan- 
tization of the transverse motion of the electron in the 
film leads to an increase of T c with decreasing film thick- 
ness, arising essentially from an enhanced effective BCS 
pairing interaction. In the previous calculations, 
the phonon modes were assumed to be the same as in 
bulk material and only the one dimensional quantum 
confinement effect of electrons were considered. Ref- 
erence 16 considered the phonon dispersion in a thin 
film undergoing substantial modification compared with 
the bulk system with quantization of phonon spectrum. 
They found the resonant shape of superconductor tran- 
sition temperature T c arises from both electronic and 
phonon confinements. Because of the quantization of 
both electron and phonon energies, the effective electron- 
electron interaction modified by the quantized phonon 
is different from the interaction arising from the bulk 
phonon. E3 They found some fine structures in the en- 
ergy gap and critical temperature by comparing with the 
bulk phonon model. [lfij| However experimentally, no evi- 
dence of such fine structure from phonon confinement was 
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found. [TlUlij Thus, we will evaluate the thermodynamic 
characteristics by the same multi-band BCS theory, but 
neglect the phonon quantization effects. 

We will start our model by taking into account of a 
realistic boundary condition. The thin film is confined 
in the x direction with geometric thickness Nt. N is the 
thin film layer and t is the average inter-layer distance. 
Because of the finite potential barriers at the interface 
of thin film with Si substrates and vacuum, we will ex- 
pand the film boundaries slightly to allow certain amount 
of charge spillage at two interfaces. |18l Il9j| For simplic- 
ity, we take charge spillage distances as the same at two 
interfaces. So the physical thickness is d = Nt + 2A(j. 

The multi-band BCS Hamiltonian of the system is 
given bv[lfij| 
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where c k is the electron operator in the nth sub-band 
with spin er, £ n (fc) = e«(fc) — [i the electron energy in the 
sub-band n measured from the chemical potential and 
^fefc™ t ne attractive interaction between nth sub-band and 



mth sub-band is given bv[l8j] 
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Only if and < Ti^d- J is the interaction 

constant and L is the periodic distance in the y and z 
directions. hu>D is the Debye energy. The constant b 
comes from the integral cutoff. 

In the weak-coupling approximation, we have the su- 
perconducting energy gap for the nth sub-band given by 
the gap function using the mean field method 
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where E mk = (^ fc 

The off-diagonal terms in the summation reflect the 
possibility of transition of the electron pair from one sub- 
band into another as the result of interaction with the 
confined phonons. Integration over k gives the gap func- 
tion of the nth sub-band at T. After some derivation, we 
arrive at the gap equation at zero temperature 



A = hojjj sinh 1 (- 
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The critical transition temperature is 
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where we assume the integral range of k' in Eq. (3) are 
the same for each subband. K is a special interaction 
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FIG. 1: (a) The critical temperature and (b) the supercon- 
ducting energy gap of the nth occupied subband versus super- 
conducting film layers N, compared with bulk values (straight 
line). Here we take Fermi wave length fc/ = 1.57A , interac- 
tion constant k = 1.5A , Debye energy huio = 0.01e«, with 
other parameters as b — 1A, t — 2.8A and Ao = 0.7A. All 
the following graphs are with the same parameters. In our 
models, the energy gap is independent of each subband. 



constant given by K — nh 2 /mJ and 7 is the highest 
occupied sub-band, given by the integer part of dkf/ir. 
The Fermi level only shifts slightly within 1% of its bulk 
value at very thin thickness. yjj So in our numerical cal- 
culation, we choose the Fermi wave vector as a constant, 
independent of film thickness. These results are similar 
as in Ref. 18. The bulk values of energy gap and transi- 
tion temperature are given as:|l8| 
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In Fig. 1, we compare the results of energy gap at 
T = OK and the superconducting transition temperature 
T c with those of bulk values. We could see that the gap 
and T c all oscillate with a period of 5 atomic layers. This 
can be understood by considering that the ratio of half 
Fermi wave length versus atomic layer distance is roughly 
2/3. If the layer thickness changes continuously, the os- 
cillation period is known to be A//2. The modlulating 
of these two quantities will lead the oscillation pattern 
differ by 5ML. We also notice that the oscillation ampli- 
tude will decay with increasing of the film thickness. But 
all these oscillation values are smaller than bulk values. 
This is in accord with several experimental results. [ill IT3 | 
And the most interesting thing is that the energy gap A„ 
oscillate the same way as T c , which leads to the universal 
constant A„/fcsT c = 1.76. This constant is the same as 
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the bulk values. Similar conclusions have been reached 
by experiments in ultra thin Pb films .[l4| but a higher 
universal constant of 4.4 is obtained. Because Pb has a 
very strong electron phonon coupling, the weak coupling 
BCS theory used here is not quite applicable. 

The Gibbs free energy, the thermodynamic critical field 
H c , the entropy and the specific heat per volume in a 
superconducting film are given, respectively by 
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FIG. 2: The thermodynamic critical field H c in the unit of 
bulk value (normalized H c ) as a function of film layers N in the 
unit of bulk values at the reduced temperature T/T c — 0.1. 
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where f km = [exp((3E km ) + l]" 1 and (3 = \/k B T. 

The thermodynamics critical field H c is also shown in 
Fig. 2. We also find the similar oscillation of H c with the 
the film layers N, although no experimental are available 
yet. The layer dependent H c shows an oscillation ampli- 
tude of 22% of the bulk situation at 5 to 8ML and decay 
to the amplitude of 5% of the bulk H c at thicker films. 
However, QSE on the perpendicular upper critical field 



FIG. 3: (a) The variation of specific heat with reduced tem- 
perature T/T c for different film thickness. The normalized 
specific heat versus film layers N at superconducting state 
(T/T c — 0.9) and at normal state (T=8k) is shown respec- 
tively in (b) and (c). 



H C 2 in the ultra-thin lead film has been reported whose 
oscillation is out of phase with the H c oscillation. |20| 

The specific heat is also calculated with Eq. (11) at 
both superconductor and normal state. In Fig. 3(a), we 
show the specific heat as a function of reduced tempera- 
ture A (A = /T c ) at different films layers. We could see 
the jump of specific heat at critical temperature. Mean- 
while, we also notice that although the specific heat both 
at normal and superconducting state oscillate with film 
thickness, the jump of specific heat at T c divided by the 
normal state specific heat [(C s — C„)/C„]t c is another 
universal constant 1.42, and this constant is the same 
as in bulk situation. In Fig. 3(b) and 3(c), we can see 
that the specific heat at supercondcting state varies in 
the same way as T c does. But its value could be slightly 
higher than the bulk case. For the specific heat in a nor- 
mal state, it oscillate unsymmetrically around its bulk 
value. This oscillation behavior is different from others. 
Hopefully, all these results could be verified in future ex- 
periments. 

The ratio of superconductor-normal tunneling differ- 
ential conductance Gsn to normal tunneling differential 
conductance Gnn is 



Gsn 
Gnn 
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The tunneling differential conductance versus bias volt- 
age calculated at different film layers and different tem- 
perature from Eq. (12) are also shown in Fig. 4. The 
peak of the conductance is almost at eV = A n . Because 
in our model all the energy gaps of each sub-band are as- 
sumed to be the same, we could only find one peak in the 
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FIG. 4: ^ SJV as the function of applied voltage V (a) at 
different film thickness with T = IK and (b) for different 
reduced temperature A with film layers N=7. 

tunneling conductance. We also find that the position of 
the conductance peak oscillates with thickness. The rea- 
son for this kind of oscillation is due to the oscillation 
of energy gap A„ at different thicknesses. In Fig. 4(b), 
we can see that the peak is depressed at high tempera- 
ture, because energy gap decreases with increasing tem- 
perature. Recently, Daejin et al. have already extracted 
energy gap from the measured conductance spectra to 
determine the transition temperature. [Tij 

In this paper, we apply multi-band BCS theory to cal- 
culate thermodynamic quantities and find their oscilla- 
tions with the film thickness, including the thermody- 
namic critical field H c , the specific heat at normal and su- 
perconducting state and the position of differential con- 
ductance peak. These oscillations are the manifestations 
of the QSE. But two universal constants A n /ksT c = 1.76 
and \{C S — C„)/C„]t c = 1-42 is independent of the film 
thickness, and their values are the same as in bulk situ- 
ation. We hope that these interesting oscillations could 
all be observed experimentally in the future. 
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